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1. Introduction 

Leptonic decays of charmed mesons were not expected to be a channel where new physics 
might be found. However, new, precise experimental results by CLEO show unexpectedly high 
rates in the decays Ds ^ TV compared to estimates from decay constants in the quenched ap- 
proximation. The HPQCD collaboration found the effect of (rooted) dynamical staggered quarks to 
be significantly smaller than the difference between experiment and the quenched calculations ' . 
Is this evidence for new physics [^], is it a statistical fluctuation or underestimate of systematics 
in the experiment, or is it a systematic effect unaccounted for by the errors quoted in We are 
aiming at a precise calculation of fjj^ as well as other observables such as the charm quark mass, 
using the Nf = 2 CLS configurations which reach small lattice spacings, where the charm quark 
mass in lattice units is really small. Here we describe first, encouraging, steps. In particular, we 
find small lattice spacing effects for 0(a) improved Wilson quarks. 

2. The CLS coordinated lattice simulations effort 

Coordinated Lattice Simulations (CLS) is a community effort to bring together the human and 
computer resources of several teams in Europe interested in lattice QCD. CLS member teams are 
located at CERN, in Germany (Berlin, DESY/Zeuthen, Mainz), Italy (Rome) and Spain (Madrid, 
Valencia). All CLS simulations use M. Liischer's implementation of the DD-HMC algorithm [Q] 
to efficiently simulate Nf = 2 Wilson QCD with non-perturbative 0(a) improvement on a variety 
of computer architectures ranging from PC clusters to the BlueGene/P at NIC/Forschungszentrum 
Jiihch. 

Table |l| shows the existing CLS ensembles. For this initial study of heavy quark physics on 
these ensembles we will use the D2, E6 and Q4 ensembles in order to get an idea of the size of the 
sea quark mass and lattice spacing effects. 

3. Setting the scale 

A final determination of the lattice scale, e.g. via niQ, is not yet available for the CLS ensem- 
bles. In ^ Del Debbio et al. determined the scale on the coarsest (j3 = 5.3) CLS ensembles to 
be a = 0.0784(10) fm via a combination of itik and niK*- We use this as our value for a on the 
D2 and E6 ensembles, and run it to jS = 5.7 for the Q4 ensemble by means of the scale L* defined 
in via g^{L*) = 5.5 in the Schrodinger functional scheme. Specifically, we use the linear fit 
log(LVa) =2.3338 + 1.4025(j8 -5.5) ±0.02. 

Since the uncertainty about the scale is an important source of error at j8 > 5.3, and the some- 
what unphysical determination of the scale may be considered as a source of an unquantifiable 
systematic error even at j8 = 5.3, a more accurate determination of the scale is certainly a priority 
in order to make accurate predictions. 

' Note that a relatively new discretization for the charm quarks is used. 



2 



Ds physics from fine lattices 



Georg von Hippel 



4. Measurements and Analysis 

We use 6 time-localized U{1) noise sources per configuration to measure the con^elators Caa, 
Cap, CpA and Cpp, where 

Cxy(xo) = -a'Y,{Xn{x)Y2im , (4.1) 

X 

Pij = qijsqj and Ajj = qiJoYsqj, on 61 configurations of D2, 28 configurations of E6, and 31 con- 
figurations of Q4, performing a fully correlated error analysis using the Jackknife method in each 
case. 

As in we define the effective mass Meff(;co) via 

g(Meff(xo),xo - a) C{xo - a) 



g{Mes{xo),xo) C{xo) 
where g{M,x) = e^^-" + e^'^^'^^^K Effective matrix elements are defined as e.g. 



(4.2) 



V ^(Meff(xo),Xo) 

We also define the PCAC quark mass as 

. , ^ .X ^{do + d*)CpA{xo)+CAadod*Cpp{xo) 
[m, + mc) = m(xo) = Cpp{xq) ' 

which needs to be renormalized and 0(a) improved as m{ji.) = Z^Zp '(/x)(l + \{bA —bp){mqa))m. 
In terms of these quantities, the (renormalized) pseudoscalar decay constant is defined as Fps = 
ZA{l + \bA{m,a))'^. 

We use non-perturbative renormalization wherever possible, in particular for ca [§], Z^ 
Zp [[lO|], and bA — bp [|ll|]. Perturbative (one-loop) renormalization is used only for bA, where no 
non-perturbative results are available. We translate the PCAC masses into the RGI masses through 
non-perturbative running in the Schrodinger functional scheme as in [[l0|]. 



Id 


Size 


a [fm] 


K 


MD T 


Id 


Size a [fm] 


K 


MD T 


Dl 


48 X 24^ 


0.08 


0.13550 


2575 


Ml 


64 X 32^ 0.06 


0.13620 


4055 


D2 






0.13590 


2565 


M2 




0.13630 


3772 


D3 






0.13610 


2520 


M3 




0.13640 


2980 


D4 






0.13620 


2505 


M4 




0.13650 


3790 


D5 






0.13625 


2510 


M5 




0.13660 


2570 


El 


64 X 32^ 


0.08 


0.13550 


2672 


PI 


96 X 48^ 0.04 


0.13620 


1702 


E2 






0.13590 


2512 


P2 




0.13630 


started 


E3 






0.13605 


2512 


P3 




0.13640 


started 


E4 






0.13610 


2497 


Q4 


128 X 64^ 0.04 


0.13640 


1450+ 


E5 






0.13625 


2656 


Q5 




0.13650 


1120+ 


E6 






0.13635 


4960 


Q6 




0.136575 


started 



Table 1: The existing and running CLS ensembles; ensembles used in this study are shown in boldface. The 
molecular dynamics time T is given in MD units after thermaUsation; trajectory length is typically T = 0.5. 
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Figure 1: The quark mass Mc (divided by 10 for scale) and fo^ as a function of (mj^ — \m\) on the E6 
ensemble, together with the Unearly extrapolated values at the physical point. 

5. Preliminary results 



5.1 Ds masses and Mc 

Using two heavy quark masses on each ensemble, and two light/strange-quark masses on E6 
and Q4, three light/strange quark masses on D2, and one additional light quark mass on E6, we 
measure the masses of all possible mass combinations of pseudoscalar mesons, as well as the 
corresponding PCAC quark masses. From the latter, we extract the RGI mass Mc as described in 
©HI- 

To extract a physical value for Mc from these data, we first interpolate Mc linearly as a function 
of Md^ to get Mc as a function of the light and strange quark masses. This we treat as a function of 
{m\ — \m\), a ;^PT-inspired proxy of the strange quark mass, and extrapolate to the physical point 
(m| - \^ml) = 0.236 (GeV)^ as shown in fig. |. 

Our results are Mc = 1694(3)(34) MeV, 1767(15)(35) MeV and 1666(1)(33) MeV on the D2, 
E6 and Q4 ensembles, respectively. We note that the lattice spacing dependence (the 1% difference 
between the D2 and Q4 ensembles) is small, but the sea quark mass dependence (the 4% difference 
between the E6 and D2 ensembles) is noticeable. 



5.2 Lattice spacing effects in fr,^ 

To get an estimate of what the lattice spacing effects on the decay constant of the Dg are likely 



to be, we define an (unphysical) reference point at nf^ = nf^ = 618 MeV, mo = 1968 MeV, to 
compare results obtained at different lattice spacings. This point is realized directly on the D2 
ensemble, where L*/^^ = 0.541(16)(11), and L*/ff^ = 0.805(12)(16). On the Q4 ensemble, we 
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Figure 2: Summary plot showing the dependence of our results on the pion (sea quark) mass, as well as the 
HPQCD result and the experimental value for comparison. 

need to interpolate in to obtain L*ff = 0.578(9)(12)(6) and L*/]^^ = 0.797(9)(16)(9). The 
errors quoted are from statistics, L* scale setting and interpolation, respectively. We find that lattice 
spacing effects are about 7(5)% in f^^, but small in fj^^. 



5.3 towards the physical point 

To approach the physical point, we take our lightest pion mass, which is nij^ = 234(10) (3) 
MeV on the E6 ensemble. 

As for the quark mass, we interpolate linearly to mjj^ = 1968 MeV at fixed {mj^ — \^m\). A 
plot of fo, as a function of (m\ — jm\) is shown in fig. |l]. Extrapolating to the physical point 
(m| - \ml) = 0.236 (GeV)2, we obtain our preliminary estimate of = 257(3)(3)(5)(?), where 
the question mark denotes unknown systematic errors, including those coming from the overall 
scale ambiguity and the quenching of the strange and charm quarks. 

We summarize our preliminary findings for /o^ in fig. ^ which illustrates that cutoff effects 
are small (at least for heavy pions), and that the light-quark mass dependence is also small (at least 
on the coarser lattice). The chiral and continuum extrapolations therefore seem to be well possible. 



6. Summary 

The CLS effort is now simulating very large and fine Nf = 2 lattices, and lattice spacings as 
small as a = 0.04 fm have become accessible, making fully relativistic charm quarks feasible. As 
simulations are progressing, lighter sea quarks are also being simulated. 
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Our preliminary study of the Dg system indicates that cutoff effects are small and under control, 
but a more precise scale determination is a priority in order to eliminate an important source of 
systematic error. 

With better statistics and more sea and valence quark masses to come, we expect to be able to 
perform an accurate determination of fjj^ in the near future. 
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